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Thermally activated escape over a potential barrier in the 
presence of periodic driving is considered. By means of novel 
time-dependent path-integral methods we derive asymptoti- 
cally exact weak-noise expressions for both the instantaneous 
and the time-averaged escape rate. The agreement with accu- 
rate numerical results is excellent over a wide range of driving 
strengths and driving frequencies. 

PACS numbers: 05.40.-a, 82.20.Mj, 82.20.Pm 
The problem of noise driven escape over a potential 
barrier is ubiquitous in natural sciences Q|. Typically, 
the noise is weak and the escape time is governed by 
an exponentially leading Arrhenius factor. This scheme, 
however, meets formidable difficulties in far from equilib- 
rium systems due to the extremely complicated interplay 
between global properties of the metastable potential and 
the noise Prominent examples are systems driven 

by time-periodic forces exemplified by strong laser 
driven semiconductor heterostructures, stochastic reso- 
nance j5j, directed transport in rocked Brownian motors 
[3|, or periodically driven "resonant activation" processes 
[ 7 1 like AC driven biochemical reactions in protein mem- 
branes. Despite its experimental importance, the theory 
of oscillating barrier crossing is still in its infancy. Pre- 
vious attempts have been restricted to weak (linear re- 
sponse) , slow (adiabatic regime) , or fast (sudden regime) 
driving In this Letter we address the most chal- 

lenging regime of strong and moderately fast driving by 
means of path-integral methods. In fact, our approach 
becomes asymptotically exact for any finite amplitude 
and period of the driving as the noise strength tends 
to zero, and comprises a conceptionally new, systematic 
treatment of the rate prefactor multiplying the exponen- 
tially leading Arrhenius factor. Closest in spirit is the 
recent work [E|, which is restricted, however, to the lin- 
ear response regime for the exponentially leading part 
and treats the prefactor with a matching procedure, in- 
volving the barrier region only. Our analytical theory 
is tested for a sinusoidally rocked metastable potential 
against very precise numerical results. Conceptionally, 
our approach should be of considerable interest for many 
related problems: generalizations for higher dimensional 
systems and for non-periodic driving forces. 

Model — We consider the overdamped escape dynamics 
of a Brownian particle x(t) in properly scaled units 



x{t) = F{x(t),t) + V2D 



(1) 



with unbiased ^-correlated Gaussian noise (thermal 
fluctuations) of strength D. The force-field F(x,i) is as- 
sumed to derive from a metastable potential with a well 
at x s and a barrier at x u > x s , subject to periodic mod- 



ulations with period T. For D — 0, the deterministic 
dynamics (Q) is furthermore assumed to exhibit a stable 
periodic orbit (attractor) x s (t) and an unstable periodic 
orbit (basin boundary) x u (t) > x s (t). 

For weak noise D, there is a small probability that 
a particle obeying (]!]) escapes from the basin of attrac- 
tion A(t) := (—oo, x u (t)\ of x s (t) and disappears towards 
infinity. For an ensemble of particles with probability 
density p{x, t), the population Pj\,(t) within the basin of 

attraction is J*" p(x,t) dx and the instantaneous rate 

of escape T(t) equals —Pj^(t)/P^(t). Apart from tran- 
sients at early times, this rate T(t) is independent of 
the initial conditions at time to. Without loss of gen- 
erality we can thus focus on x(to) = x s (to). Small D 
implies rare escape events, i.e. the deviation of fU(t) 
from its initial value P^(to) = 1 is negligible. Exploit- 
ing x u (t) = F(x u (t),t) and the Fokker-Planck-equation 
dtp = d x {—F(x, t) + Dd x } p governing p — p(x, t) we find 
for the instantaneous rate 



r(t) = -Dd xP (x = x u (t),t) 



(2) 



Path-integral approach — With the choice p(x, to) — 
S(x — x s (to)) we obtain for the conditional probability 
p(x,t) the path integral representation H 



p(x f ,t f ) = / Vx(t) e - s ^/ D 



where the "action" is given by 

S[x(t)] := / dt [x(t) - F(x(t),t)} 2 /4 



(3) 



(4) 



and where x(to) = x s (to) and x(tt) = Xf are the "initial" 
and "final" conditions for the paths x(t). For weak noise, 
the integral (|^) is dominated by a set of paths x* k {t), 
corresponding to minima of the action (^) (distinguished 
by the label k) . These satisfy an Euler-Lagrange equation 
equivalent to the following Hamiltonian dynamics 



pt{t) = -pt{t)F'{x* k {t),t) 
xUt)=2pUt)+F(x* k (t),t) 



(5) 
(6) 



with F'(x,t) := d x F(x,t). For well separated paths 
x k (t), a functional saddle point approximation in (^) 
yields 



-S[xl(t)]/D 

- [AnDQKtf)} 1 / 2 



[1 + 0(D)] , (7) 



where the quantity Q* k {t) relates to the determinant of 
fluctuations around x j (t) . Following the reasoning in || , 
we find for our case |10| that Q* k {t) obeys the relation 



1 



Q* k (t)/2-d[Qt(t)F'(xt(t),t)]/dt 
+ Ql(t)p* k (t)F"(x* k (t),t) = 







(8) 



with initial conditions Q k {to) — and Ql(to) = 1- Ex- 
ploiting that the derivative of the action at its cnchpoint 
equals the momentum p k (tf), we can infer from 
first main result, namely 



(9) 



where the boundary conditions x k (tq) = x s (to) and 
xt(tf) — x u (tf) are understood in ([s],^]). In view of 
(Q), the instantaneous rate (0) has the suggestive form 
of probability at the separatrix times "velocity" . 

Closer inspection of (||||) reveals the following generic 
features of each path x k (t) which notably contributes to 
the rate ([)]), see fig. 1: Starting at x k (t$) = x s (tQ), it 
continues to follow rather closely the stable periodic or- 
bit x s (t) for some time. At a certain moment, it crosses 
over into the vicinity of the unstable periodic orbit x u (t) 
and remains there for the rest of its time, ending at 
xl(tf) = x u (tf). Without loss of generality, we can sort 
the paths x* k (t) by the time they spend near the unstable 
periodic orbit, such that XqW * s that path which crosses 
over from x s (t) to x u (t) at the "latest possible moment" . 
Apart from a time shift, each path x k (t) then closely re- 
sembles the same "master path" x*(t) (see fig. 1). This 
path x*(t) is defined as an absolute minimum of the ac- 
tion (|J) in the limit to — > — oo, tj — > oo, and is fixed 
uniquely by demanding that x* (t + kT) is the "master 
path" associated with x k {t). 

The basic qualitative features of each minimizing path 
x k (t) are thus quite similar to the well-known barrier- 
crossing problem in a static potential jnj. However, in 
the limit to — > — oo, tj -> oo we have, in contrast to 
this latter situation, not a continuous symmetry (Gold- 
stone mode) , but a discrete degeneracy of the minimizing 
paths. As a consequence, in our case the minimizing 
paths xl(t) remain well separated and thus the rate for- 
mula (|9|) is valid for any (arbitrary but fixed) finite values 
of the driving amplitude and period, provided the noise 
strength D is sufficiently small. On the other hand, for a 
given D we have to exclude extremely small amplitudes 
and extremely long or short periods since this would lead 
effectively back to the static case. 

As long as the "master path" x* (t) remains sufficiently 
close to the stable periodic orbit x s (t) , say for t < t s , the 
force-field is well approximated by 

F(x,t) = F(x s (t),t) + (x-x s (t))F'(x s (t),t) . (10) 

An analogous approximation for F(x, t) is valid while 
x*(t) remains in a sufficiently small neighborhood of 
x u (t), say for t > t u . The corresponding local solutions 
of the Hamilton equations (pLp[) can then be written as 



(11) 
(12) 



p*(t)=p*(t s , u )e- A '^ t > t ^ 
x*(t)=x s , u (t)±p*(t)I s , u (t) . 



Here 's,u' means that the index is either 's' or 'u' and 
the upper and lower signs in ( |l2| ) refer to 's' and 'u\ 
respectively. Further, we have introduced 



As,it(^, t s _ u ) 



F'(x SiU (i),t)di 



2 / 

=poo 



2A„ 



,(t,l) £ 



(13) 
(14) 



Similarly, the local solutions for the prefactor in (ph can 
be written as 



Q*(t < t s ) = I s (t)/2 

Q*(t > t u ) = Cl / P *(tf 



C2l U (t) ■ 



(15) 
(16) 



The parameters p*{t s .u) in (jllj) and Ci,2 in ( |16| ) cannot 
be fixed within such a local analysis around x StU (t), they 
require the global solution of (j|J^,^). We furthermore 
observe that due to the time-periodicity of F(x,t) and 
x s ,uit), the quantities 



K,u :=A StU (t + T,t)/T 



(17) 



are indeed t- independent. The stability/instability of the 
periodic orbits x SlU (t) implies X s < and X u > 0. It 
follows that I s . u (t) from (|14| ) are finite, T-periodic func- 
tions. 

The expressions for x* k {t), p* k {t) 1 and Q* k {t) are some- 
what more complicated than in ( pd[]l^ ) but since x k (t) 
is well approximated by x*(t + kT), the same follows 
for p* k {t) and Q* k {t). Closer inspection shows [pi that 
in (ph the pre-exponential factors p k (tf) and Q* k \tf) can 
be approximated by p*(tf + kT) and Q*Uf + kT) with- 
out further increasing the error 0{D) in (|5j). Within this 
same accuracy, the exponential in (|9|) requires - due to 
the small denominator D - a somewhat more elaborate 
approximation, yielding 



S[xt(t)} = S[x*(t)} 



p*(t) 2 dt 



(18) 



t,+kT 



Rate formula — By introducing these approximations 
into (£)]), exploiting dTl]-|l6|), and dropping the index '/' 
oft/, we obtain |1(| as the central result of this work the 
instantaneous rate 



T(t) = VDae~ 



S[x*(t)]/D 



n(t,D) [l + E(D)} 



a := [4ttT 2 lim p*{t) 2 Q*{t)]^ l/2 

t— *OG 

Mt) 2 



n(t,D) := T 
h{t) 



E 

k— — oo 



k (tyi u (t)/2D 



D 



lim p*(t) e A " ( *' 4) 



(19) 
(20) 

(21) 

(22) 



The relative error E(D) is found to be of order O(D) 
if F"(x u (t),t) ee 0, and O(Z) 1 / 2 ) otherwise. By use of 
( |l3| - |l7| ) one finds that the average of ([H]) over a single 
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time-period T equals 1. For the time- averaged rate T we 
thus obtain 



T = VDae 



,-S[x*{t)]/D 



[l + E{D)] 



(23) 



It consists of an Arrhenius-type exponentially leading 
part and, in contrast to equilibrium rates Q , a non-trivial 
pre-exponential D-dependence. 

An archetype example — In general, the explicit quan- 
titative evaluation of S[x*(t)], a, and n(t,D) in ( |19|]23| ) 
is not possible in closed analytical form. An exception 
is the piecewise linear force-field with additive sinusoidal 
driving 



F(x <0,t) = X s (x — x a )+A sin(fit) 
F(x > 0,i) = X u (x - x u ) + A sin(fH) , 



(24) 



corresponding to a periodically rocked piecewise 
parabolic potential, with parameters x s < 0, x u > 0, 
A s < 0, A„ > 0, respecting X s x s — X u x u (continuity at 
x = 0). To simplify the analytical calculations we further 
restrict ourselves to the case that the "master path" x* (t) 
crosses the matching point x — in (|24|) only once p2[ , 
say at t = t c . The periodic orbits x s>u (t) then assume 
the simple form 



.(*) 



A [X StU sin(f2f) + ficos(fii)] 



\2 



n 2 



(25) 



Moreover, Eqs.( |ll| , [12D are now valid with index 's' for all 
t < t c and with 'u' for all t > t c . By matching these 
solutions at t = t c the global parameters p*(t s , u ) in (JlT 
are fixed and one obtains 



i ( X U X S — 

— arctan — — — 

o \n(x u + Xs) 



(26) 



To ensure that x* (t) is a minimum of the action in 
one has 



(A« + X s ) Aftcos(flt c ) < 



(27) 



In the same manner, the prefactor ( |l6| ) has to be matched 
at t = t c - While Q*(t) is still continuous, Q*(t) develops 
a jump at t = t c which can be determined from Upon 
collecting everything, the final result reads 



S[x*(t)] = AV 

'\A\(n 2 + 



i - 



\AX U X S 



(RAV) 1 / 2 
X U X S ) + (RAV)V 2 



l6ir 3 \A\S[x*{t)} 



1/2 



Pk(t) 



-\ u (fcT +t-t c ) 



\AX U X S 



HV 2 



(28) 

(29) 
(30) 



where H := (A* + Q 2 )(X 2 S + fi 2 ), R := 2H/(X- 1 - Xj 1 ), 
and Ay := (A u a; 2 — X s x 2 )/2 is the potential barrier corre- 
sponding to the undrivcn (A = 0) force-field (PJ). With 



T = 2n/n, I u {t) = 1/A„, and t c from (g§g7|), the rate 
(|l^,|3|) is thus determined completely. 

Comparison — These analytical predictions for the in- 
stantaneous rate ( |l9| ) are compared in fig. 2 for a rep- 
resentative set of parameter values with very accurate 
numerical results. The agreement indeed improves with 
decreasing noise-strength D. While the absolute values 
of T(t) and the location of the extrema strongly depend 
on D, the overall shape changes very little and does not 
develop singularities as D — > 0. The corresponding time- 
averaged rates ( p3| ) are depicted in fig. 3, exhibiting ex- 
cellent agreement between theory and numerics even for 
relatively large D. The inset of fig. 3 confirms our pre- 
diction that the relative error E(D) in ( p3| ) decreases 
asymptotically like D. Finally, fig. 4. illustrates the de- 
pendence of the averaged rate T upon the amplitude A of 
the periodic driving force. As expected, our theoretical 
prediction compares very well with the (numerically) ex- 
act rate, except for very small driving amplitudes A (see 
the discussion above eq. ([n])). The approximation from 
H is complementary to ours in that it is very accurate 
for small A but develops considerable deviations with in- 
creasing A. Those approximations have been omitted in 
figs. 2, 3 since they are not valid in this parameter regime 
and indeed are way off. 
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FIG. 1. Solid: The paths x%(t), k = 0, . . . , 3, which mini- 
mize the action (^) with boundary conditions x%(to) = x a (to) 
and x* k {tf) = x u (tf) for the example ( pi[ ) with x a = X 3 = — 1, 
x u = X u — Q = 1, A — 0.5 (dimensionless units). Dashed: 
the associated "master paths" x*(t + kT). Dotted: stable 
and unstable periodic orbits x a (t) and x u (t). In this plot, 
tf — to has been chosen rather small. As tf — to increases, 
more and more intermediate paths x%(i) appear which better 
and better agree with x*(t + kT). 
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FIG. 3. Arrhenius plot of the time-averaged rate F. Pa- 
rameters are like in fig. 2. Solid: analytical result ( ]2^ , p8|j29[ ) . 
Crosses: precise numerical results. Inset: relative difference 
between analytical (F) and numerical (F nU m) rate. 




FIG. 4. Time-averaged rate F vs. driving amplitude A 
for D — 0.05. Parameters are like in fig. 2. Crosses: precise 
numerical results. Solid: analytical result ( p3|j2^ , p9[ ) . Dotted: 
theoretical approximation from Ref. [4]. The small- A regime 
is magnified in the inset. 



FIG. 2. Instantaneous rate V(t) for the force field (|2J) 
with x a = X a — — 1, x u = X u = Q — 1, A = 0.5, cor- 
responding to a static (A — 0) po tential barrier AV = 1. 
Solid: analytical result ( |l9|j2lH26| - |3C| ). Dashed: high- precision 
numerical results by evolving the Fokker-Planck-equation for 
p(x,t) until transients have died out and then evaluating (bl). 
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